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Abstract: Birds and other animals are frequently killed by cars, causing the death of many million individuals per 
year. Why some species are killed more often than others has never been investigated. In this work hypothesized that risk 
taking behavior may affect the probability of certain kinds of individuals being killed disproportionately often. 
Furthermore, behavior of individuals on roads, abundance, habitat preferences, breeding sociality, and health status may 
all potentially affect the risk of being killed on roads. We used information on the abundance of road kills and the 
abundance in the surrounding environment of 50 species of birds obtained during regular censuses in 2001—2006 in a rural 
site in Denmark to test these predictions. The frequency of road kills increased linearly with abundance, while the 
proportion of individuals sitting on the road or flying low across the road only explained little additional variation in 
frequency of road casualties. After having accounted for abundance, we found that species with a short flight distance and 
hence taking greater risks when approached by a potential cause of danger were killed disproportionately often. In 
addition, solitary species, species with a high prevalence of Plasmodium infection, and species with a large bursa of 
Fabricius for their body size had a high susceptibility to being killed by cars. These findings suggest that a range of 
different factors indicative of risk-taking behavior, visual acuity and health status cause certain bird species to be 


susceptible to casualties due to cars. 
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The impact of humans on animals has intensified human populations and the resulting exploitation of 
dramatically during the last decades due to increasing natural resources across the globe (Marzluff et al, 2001). 
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Urbanization has a number of consequences for wild 
animals because of the longer growing season associated 
with altered micro-climatic conditions, a change in 
abundance of habitats, an increase in food abundance and 
hence population density of specific animals, and 
proximity of humans that may be problematic for species 
with low thresholds of fear (Batten, 1973; Klausnitzer, 
1989; Gliwicz et al, 1994; Stephan, 1999; Shochat et al, 
2006; Møller, 2008a). A major impact of urbanization is 
the construction of roads that disrupt animal movement, 
but also affect everyday life of animals due to the 
presence of cars and their effects. For example, recent 
research has shown strong effects of noise from traffic on 
song learning and change in composition and carrier 
frequency of songs (e. g. Reindt, 2003; Warren et al, 
2006). Second, roads can have dramatic effects on the 
distribution of animals (Reijnen & Foppen, 1991; 
Underhill & Angold, 2000; Forman & Alexander, 1998). 
A third consequence is the large number of traffic kills 
due to cars (Blimel & Bliimel, 1980; Bruun-Schmidt, 
1994; Dunthorn & Errington, 1964; Erritzøe et al, 2003; 
Finnis, 1960; Göransson et al, 1978; Haas, 1964; 
Heinrich, 1978; Hodson, 1960; Hodson & Snow, 1965; 
Nankinov & Todorov, 1983; Smettan, 1988; Straube, 
1988; Vignes, 1984; Wäscher et al, 1988). For example, 
the total number of birds killed by traffic in Denmark 
alone (a small country of only 43,000 km’) was 
estimated to be 3.50 million in the 1960’s and 3.27 
million in 1980 (Hansen, 1969, 1982). Therefore, the 
global estimates must at least be hundreds of millions per 
year. 

While the effect of cars on avian mortality is well 
described, there is to the best of our knowledge no 
information about the causes of interspecific variation in 
such mortality. Superficially, one might assume that all 
species are equally susceptible to traffic mortality 
because the recent invention of cars has not allowed 
adaptation to this cause of mortality. However, there 
might be behavioral or physiological traits that 


predispose particular species to elevated risk of mortality. 


Are road-kills just reflecting a random sample of all 
species according to their abundance and the probability 
that individuals sit on the road or fly close to the ground, 
thereby exposing them to traffic accidents, or do specific 
characteristics of certain species render them more 
susceptible to being killed by cars? The objectives of this 
study were to address to which extent the abundance of 
road-kills of different species could be explained by 


ecological factors. To this end we estimated both the 
frequency of road-kills of different species, but also their 
expected frequency based on standardized birds censuses. 

The aims of this study were to assess risk factors 
associated with road-kills of birds. As a null hypothesis 
we assessed the importance of abundance as a predictor 
of the frequency of road-kills of different species, 
because abundant species by necessity are more likely to 
end up being killed by cars. Next, we calculated an index 
of susceptibility to mortality due to cars by quantifying 
residual variation for each species. First, we 
hypothesized that species living in open habitats would 
have a better ability to escape an approaching car 
because of a better view of the approaching vehicle. 
Second, Blumstein (2006) has in a study of the 
evolutionary ecology of fear suggested that flight 
initiation distance when approached by humans provides 
a reliable metric of wariness. Flight distance to an 
approaching human can be considered a metric of the 
risk that individuals of a given species are willing to take 
when encountering a potential predator (Blumstein, 
2006). Importantly, Moller et al (2008) found that bird 
species with relatively short flight distances suffered 
from higher susceptibility to predation by sparrowhawks 
Accipiter nisus, and susceptibility of birds to cat 
predation is also related to flight distance (Meller et al, 
2010). Therefore, species with short flight distances 
should be particularly susceptible to traffic. Species 
feeding on fast-moving prey may be better at motion 
detection (Fleishman et al, 1995) and have more acute 
visual systems (Garamszegi et al, 2002), making it likely 
that species eating moving prey would be more flighty 
than herbivores. Hence, species eating moving prey 
should be less susceptible to traffic than herbivores. 
Sociality may involve allocation of more time to monitor 
conspecifics (Roberts, 1988), making social species more 
vigilant to an approaching human being (Blumstein, 
2006), but incidentally also to approaching cars. Third, 
we predicted that a relatively short duration of the 
nestling period for a given body size would imply that 
fledglings had not fully developed when leaving the nest 
(Starck & Ricklefs, 1998), resulting in an increased risk 
of mortality due to cars in such species. Fourth, we 
predicted that malaria and other blood parasites would 
affect metabolic efficiency of hosts (Atkinson & Van 
Riper, 1991; Chen et al, 2001), thereby reducing working 
ability and hence escape ability when a car is rapidly 
approaching. Thus, we predicted that bird species with a 
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high prevalence of blood parasites, especially intra- 
cellular blood parasites that cause erythrocytes to burst 
and hence give rise to anemia, would be positively 
related to susceptibility to car fatalities. Finally, we 
expected that aspects of immunity in young birds would 
predict susceptibility to car accidents because such 
measures of immunity would have evolved in response 
to parasite-mediated selection, causing species with 
virulent parasites to have particularly large immune 
defense organs for their body size. Thus, we predicted a 
positive relationship between susceptibility to car 
accidents and relative size of the bursa of Fabricius, 
which is an important immune defense organ involved in 
production of antibodies and responsible for 
differentiation of the repertoire of B-cells in young birds 
(Glick, 1983, 1994; Toivanen & Toivanen, 1987). Many 
studies have shown that relative size of the bursa of 
Fabricius is related to parasite impact on avian hosts (e.g., 
Garamszegi et al, 2003, 2007; Garamszegi & Moller, 
2007; Møller, 2009; Møller & Erritzoe, 1996). 

In the present study we were able to assess traffic 
casualties in relation to null expectations based on 
In addition, we 


abundance and flight behavior. 


performed post mortem autopsies of almost all 
individuals, verifying that they indeed had died from 
impact with a car. We tested our predictions concerning 
risk factors associated with traffic fatalities of birds by 
relying on an extensive data set of road-kills collected by 
JE and HE, who also conducted standardized bird 
censuses in the same sites to quantify the abundance of 
different species. Therefore, the present study assessed 
risk factors of traffic mortality in birds, because we could 
compare observed and expected abundance of road-kills 
of different species and relate such an index of 
susceptibility to flight behavior, and ecological and life 
history characteristics. Thus, we were able to assess the 
relative importance of a range of different predictors of 


susceptibility to traffic accidents. 


1 Materials and Methods 


1.1 Road kills 

JE and HE biked or walked 10,500 km on a fixed 
route during 2001—2006 around Taps (55°27’ N, 9°26’ E), 
Denmark, with an even distribution of these sampling 
events throughout the year. The habitat of this open 
farmland area is open fields, with scattered hedges, 
plantations and houses with gardens. Traffic is not 
intense because less than 70 cars were recorded on most 


censuses. Each trip covered at least 6 km, and they found 
a total of 461 road-killed birds by carefully searching the 
road and the verges during these trips. All individuals 
were aged according to standard criteria reported by 
Svensson (2006). In addition, all individuals were 
inspected post mortem, and traffic casualty rather than 
other causes of mortality was ascertained from severe 
fractures, blood-filled lungs and other signs of a forceful 
impact between a bird and a hard obstacle (a car). Thus, 
we could exclude the possibility that the birds recorded 
had died for other reasons such as collision with 
overhead wires, poisoning, predation or other causes of 
death. The species and their abundance are listed in the 
Appendix. 

JE also recorded all individuals seen during the 
regular trips, carefully recording whether an individual 
was sitting on the road or road verge, and whether an 
individual flew at a height of less or more than 2 m, as 
judged from our own height. We used this information to 
estimate the proportion of individuals of each species 
that was sitting on the road and the proportion of 
individuals that flew below a height of 2 m as indicators 
of elevated risk of impact with cars. 

A total of 25 of 249 sparrows Passer sp. (10%) 
could not be assigned to species with certainty, and we 
attributed these unidentified sparrows to house sparrows 
Passer domesticus and tree sparrows Passer montanus 
according to the relative abundance of these two species 
as road kills of actually identified individuals. However, 
the results remained unchanged if these 25 unidentified 
individuals were excluded from the analyses. 

1.2 Bird censuses 

JE conducted standardized point counts (Bibby et al, 
2005) to quantify the abundance of birds during the 
breeding season and the winter along the roads where 
road-kills were recorded. In brief, the point counts 
consisted of 5 minutes being spent per point looking and 
listening for birds at any distance from the point, while 
recording these, with 100m between individual census 
points, using a total of 20 points per census. Total 
abundance during the breeding season was estimated as 
the sum of the seven annual point count censuses made 
between 2001 and 2006 during the breeding season 
between 25 May and 1 June. Two winter censuses were 
conducted on 5 January 2005 and 5 January 2006. We 
obtained an estimate of total abundance of each species 
by adding the number of breeding season observations to 
the number of winter observations multiplied by 3.5, to 
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account for the fact that there were 3.5 times as many 
breeding season censuses as winter censuses. The bird 
count results were repeatable among censuses during the 
breeding season (one-way ANOVA based on the logio- 
transformed counts: F = 3.72, df = 33,170, P < 0.0001) 
and winter (F = 4.66, df= 7,8, P = 0.023), implying that 
the abundance of different species was consistent among 
years. The results of these censuses are reported in the 
Appendix. 

1.3 Ecological variables 

We scored the main breeding habitats as open 
grassland (a score of 0), bushes (a score of 1), or trees (a 
score of 2), relying on information in Cramp & Perrins 
(1977-1994). 

During March-August 2006-2007 APM estimated 
flight distances for birds in Ile-de-France and Bretagne, 
France and Northern Jutland, Denmark, using a standard 
technique developed by Blumstein (2006). Estimates of 
flight distances are highly consistent among countries, as 
reported by Moller (2008a, b), hence justifying the use of 
flight distance data from another population than that 
used for studying traffic mortality. In brief, when an 
individual bird had been located while resting, foraging, 
preening or singing, APM moved at a normal walking 
speed towards the individual, while recording the 
number of steps (which approximately equals the number 
of meters). The distance at which the individual took 
flight was defined as the flight distance, while in 
accordance with Blumstein (2006) we defined starting 
distance as the distance from where the observer started 
walking until the position of the bird. If the individual 
was positioned in the vegetation, the height above 
ground, when APM started walking towards it, was 
recorded to the nearest meter. While recording these 
distances, APM also recorded date and time of day, age 
and the sex of the individual if external characteristics 
allowed sexing with binoculars. Juveniles were only 
included when they were fully independent and not fed 
by adults. Flight distance was estimated as the square- 
root of the sum of the squared horizontal distance and the 
squared height (Blumstein, 2006). 

All recordings were made during the breeding 
season, when most individuals are sedentary, thus 
preventing the same individual from being recorded in 
different sites, with each site being more than 100m 
apart (equaling a territory of 3.1 ha), and for larger 
species such as herons and corvids APM used a distance 
of 500 m (equaling a territory of 78.5 ha) because of 


Vol. 32 


their larger territories. Using only breeding birds may 
cause another problem if territorial individuals are less 
likely to have long flight distances because they are 
reluctant to cross territory boundaries. However, we can 
dismiss this alternative explanation because an analysis 
of 24 species (based on 2018 observations), assessed 
during winter (December 2006 — February 2007) and 
summer revealed a strong positive relationship between 
the two sets of estimates (weighted regression: F = 41.50, 
df= 1, 3, 17 = 0.64, P< 0.0001). 

We assessed the spatial and temporal consistency in 
estimates of mean flight distances in three ways and 
found high degrees of consistency (Meller, 2008a,b). 

Breeding sociality was quantified as colonial (a 
score of 1) when individuals used territories that only 
contained a nest site, and otherwise solitary (a score of 0), 
using information in Cramp & Perrins (1977-1994). 

We extracted information on duration of the 
incubation and the nestling period from Cramp & Perrins 
(1977-1994), relying on the value based on the largest 
sample size if more than a single value was reported. If a 
range of values was reported, we used the mean of the 
extreme values. Likewise, we extracted mean body mass 
from the breeding season from Cramp & Perrins 
(1977-1994), using the largest sample if more data were 
reported, and using the mean of means for adult males 
and adult females, if masses were reported separately for 
the two sexes. 

We used information on prevalence of four genera 
of blood parasites (Leucocytozoon, Haemoproteus, 
Plasmodium, Trypanosoma), assuming that this would 
reflect prevalence in the population, as done in previous 
studies (e. g. Moller & Nielsen, 2007). We also extracted 
information on the number of individuals examined for 
each of the host species. In total, the study was based on 
examination of infection level of 16,995 individual 
juvenile and adult hosts belonging to 263 species based 
on blood smears, with a range from 2 to 1539 individuals 
each for the species [see Moller & Nielsen (2007) for a 
complete list of sources]. Most of the blood parasite 
information derived from Northern Europe, where the 
study of road-kills was also conducted, providing a high 
degree of overlap between the geographical location of 
study sites for parasites and car casualties. Scheuerlein & 
Ricklefs (2004) have shown that prevalence estimates are 
repeatable across study sites, despite considerable 
variation among sites and sampling periods. Such 
variation among sites and sampling periods would make 
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any investigation of relationships between prevalence 
and ecological variables conservative. 

We used information on the mass of the bursa of 
Fabricius from Meller et al (2005) and Garamszegi et al 
(2007). 

The entire data set is reported in the Appendix. 

1.4 Statistical methods 

We estimated a susceptibility index to quantify 
likelihood for individuals of a given species to be killed 
by traffic, as logio-transformed number of road-kills 
minus logio-transformed (number of individuals 
observed of that species divided by the total number of 
individuals of all species) multiplied by the total number 
of road-kills). To avoid problems with species having an 
abundance of road-kills or an overall abundance of zero, 
we added one to all observations. 

Nestling period, mass of bursa of Fabricius, body 
mass and flight distance were log)o-transformed, while 
the proportion of individuals sitting on the road, the 
proportion of individuals that flew at a height of less than 
2 m and prevalence of blood parasites were square-root 
arcsine-transformed before analysis. 

Most statistical approaches assume that each data 
point provides equally precise information about the 
deterministic part of total process variation, i.e. the 
standard deviation of the error term is constant over all 
values of the predictor variables (Sokal & Rohlf, 1995). 
We weighted each observation by sample size (the 
number of individuals killed by traffic) in order to use all 
data in an unbiased fashion, thereby giving each datum a 
weight that reflects its degree of precision due to 
sampling effort (Draper & Smith, 1981; Neter et al, 
1996). Comparative analyses may be confounded by 
sample size if sampling effort is important, and if sample 
size varies considerably among taxa. In order to weight 
regressions by sample size in the analysis of contrasts, 
we calculated weights for each contrast by calculating 
the mean sample size for the taxa immediately subtended 
by that node in the phylogeny (Møller & Nielsen, 2007). 

We calculated an estimate of effect size as the 
Pearson product-moment correlation, using the equations 
in Rosenthal (1991, pp. 73-74), relying on Cohen’s (1988) 
conventions where r = 0.1 equals a small effect, r = 0.3 
an intermediate effect and r = 0.5 equals a large effect). 

The three most common species accounted for 75% 
of all road-kills, implying that these species may have a 
disproportionate influence on the conclusions. We re-did 
all analyses after exclusion of these three common 


species, but all conclusions remained unchanged, 
suggesting that the three most common species did not 
bias the findings. 

1.5 Comparative analyses 

We controlled for similarity in frequency of traffic 
casualties among species due to common ancestry by 
calculating standardized independent linear contrasts 
(Felsenstein, 1985), using the computer program CAIC 
(Purvis & Rambaut, 1995). Standardization of contrast 
values was checked by examination of absolute values of 
standardized contrasts versus their standard deviations 
(Garland et al, 1992). Plotting the resulting contrasts 
against the variances of the corresponding nodes 
revealed that these transformations made the variables 
suitable for regression analyses. In cases where extreme 
residuals were recorded, we tested for the robustness of 
the conclusions by excluding contrasts with studentized 
residuals greater than 3.00 (Jones & Purvis, 1997). 
Likewise, we ranked the independent variable to test if 
the conclusions remained unchanged (Møller & Birkhead, 
1994), and in no case did this procedure give rise to 
conclusions different from those obtained with the 
contrast values. Ranking provides a very conservative 
test of a hypothesis, and robustness of findings to 
ranking of the independent variable thus suggests that 
distributions of variables are not a confounding factor 
leading to specific conclusions. 

We constructed a composite phylogeny of the 50 
Species investigated to calculate standardized linear 
contrasts, mainly based on Sibley & Ahlquist (1990), 
combined with information from more recent sources 
(Blondel et al, 1996; Barker et al, 2001, 2004; Crochet et 
al, 2000; Helbig & Seibold, 1999)(Fig. 1). Because this 
composite phylogeny was derived from different studies 
using different phylogenetic methods and different 
molecular markers, consistent estimates of branch 
lengths were unavailable. Therefore, we considered 
branch lengths to be equal in the analyses (this is 
equivalent to an assumption of a punctuated model of 
& Rambaut, 1995). We used 
MacClade to construct the phylogeny and make 


evolution (Purvis 


inferences about the ancestral state of susceptibility to 
traffic (Maddison & Maddison, 1992). 

Regressions of standardized linear contrasts were 
forced through the origin because the comparative 
analyses assume that there has been no evolutionary 
change in a character when the predictor variable has not 
changed (Purvis & Rambaut, 1995). We used multiple 
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Fig. 1 Phylogeny of species of birds analyzed for road-kills 


Species with a relatively low (open branches) and high frequency of road-kills (dark branches), relative to expectations from bird census data, are shown. 


Hatched branches have equivocal states. See Material and methods for sources. 


regression to find the best-fit model, using the software 
JMP (Anon, 2000). The best-fit model was determined 
using Akaike’s information criterion as an estimate of 
the improvement in fit for addition of variables 
(Burnham & Anderson, 2002). There was no evidence of 
collinearity between variables (Tabacknik & Fidell, 
1996). 

Information was missing for some species for 
certain variables, causing sample sizes to differ slightly 
among analyses. 

Values reported are means (SE). 


2 Results 


The number of individuals killed varied among 50 
bird species from 0 to 199, with a mean = 9.2, median = 
2.0, SE = 4.4. The three most abundant species were 
house sparrow (199 individuals), blackbird Turdus 
merula (95) and tree sparrow (49). The total abundance 
during censuses of the 50 species varied from 0 to 261 
individuals, with a mean = 12.5, median = 2.5, SE = 5.4. 
The frequency of traffic kills increased significantly with 
abundance in a linear fashion, explaining 42% of the 
variance (Fig. 2), and the increase was significant both 
for analyses based on species-specific data and 
independent contrasts (Tab. 1). The susceptibility index 
varied from —1.235 to +0.699 with a mean of 0.005 (SE 
= 0.061), not differing significantly from zero (one- 


sample t-test, £ = 0.08, df = 49, P 0.94). The 
distribution of the relative frequency of road-kills on the 
phylogeny showed that the ancestral state among the 
species studied was a low frequency of traffic casualties 
(Fig. 1). There were many transitions between different 
states (Fig. 1), 
associated with a relatively high frequency of road-kills 


suggesting that ecological factors 


were labile. 
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Fig. 2 Frequency of traffic kills in different species in 
relation to their abundance 


The line shows the linear regression line. 


The proportion of individuals flying below 2 m 
ranged from 0 to 1 among species, with a mean value of 
0.36 (SE = 0.05). The proportion of individuals on the 
road varied from 0 to 1 among species, with a mean 
value of 0.31 (SE = 0.05). These two variables were not 
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Tab. 1 Linear regression of log-transformed number of road-kills for different species of birds in relation to 


log-transformed abundance 





Variable Sum of squares df F P Slope (SE) Effect size 
Model A (Species): 
log (Abundance) 441 1 29.76 < 0.0001 0.513 (0.094) 0.62 
Error 7.12 48 
Model B (Contrasts): 
log (Abundance) 0.63 1 18.79 < 0.0001 0.452 (0.104) 0.53 
Error 1.61 48 


Model A was based on species-specific data, while model B was based on standardized linear contrasts. Effect size was estimated as the Pearson product 


moment correlation coefficient. 


significantly correlated (Pearson r = 0.06, t = 0.44, df = 
48, P = 0.66). The proportion of individuals on the road 
explained a significant amount of variation in road-kills 
(F = 6.65, df= 1,48, P = 0.013, slope (SE) = 0.32 (0.12), 
effect size = 0.35), while that was not the case for the 
proportion of individuals flying below 2 m (F = 0.38, df 
= 1,48, P = 0.54, slope (SE) = 0.08 (0.12), effect size = 
0.09). Neither variable explained a significant amount of 
variation in susceptibility to traffic (proportion of 
individuals on road: partial F = 0.81, df= 1,47, P = 0.37, 
effect size = 0.13; proportion of individuals flying lower 
than 2 m: partial F = 0.07, df= 1,47, P = 0.80, effect size 
= 0.04). Flight distance did not significantly predict 
susceptibility to traffic in a model that included body 
mass (partial F = 1.10, df= 1,42, P = 0.30, effect size = 
0.16). 

Age did not affect susceptibility to traffic because 
the proportion of yearlings was not a significant 
predictor in a model that included body mass as an 
additional predictor variable (partial F = 0.14, df= 1,33, 
P=0.71, effect size = 0.06). 

Susceptibility to traffic increased from open habitats 
and bushes to trees in analyses of species-specific data, 
but not in analyses of contrasts (Tab. 2 models A-B). 

Breeding sociality reduced susceptibility to traffic, 
after inclusion of body mass as an additional predictor 
variable (Fig. 3A; partial F = 11.52, df = 1,46, P = 
0.0014, effect size = 0.45). In contrast, food did not 
significantly predict traffic, 
inclusion of body mass as an additional predictor 
variable (partial F = 0.00, df= 1,46, P = 0.99, effect size 
= 0.01). 

Bird species with a long nestling period for their 


susceptibility to after 


body size were less susceptible to traffic, and this effect 
was independent of the effect of habitat in analyses of 


species-specific data and contrasts (Tab. 2 Models A-B). 
In contrast, the duration of the incubation period was not 
a significant predictor (results not shown). In addition, 
susceptibility increased with body mass, although only in 
analyses of species-specific data (Tab. 2 Models A-B), 
implying that large species were more susceptible than 
small species in species-specific analyses. 

The prevalence of two genera of blood parasites, 
Leucocytozoon and Plasmodium, were both positively 
related to susceptibility to traffic, with similar effects in 
analyses of species-specific data and independent 
contrasts (Tab. 2 Models C-D; Fig. 3B). In contrast, the 
prevalence of Haemoproteus and Trypanosoma was not a 
significant predictor (results not shown for simplicity). 

The size of the bursa of Fabricius was positively 
related to susceptibility to traffic, both in analyses of 
species-specific data and independent contrasts (Tab. 2 
Models E-F; Fig. 3C). This effect was independent of 
body mass. 

Finally, we developed best-fit models that included 
all the predictor variables. An analysis of species- 
specific data revealed a model that explained 60% of the 
variance (Tab. 3 Model A). There were statistically 
significant effects of flight distance after inclusion of the 
other variables of interest (Fig. 3D; negative relationship 
between distance and susceptibility), breeding sociality 
(solitary species having greater susceptibility than 
colonial species), and mass of bursa of Fabricius (species 
with larger bursa being more susceptible) (Tab. 3 Model 
A). A model based on independent contrasts revealed 
intermediate to strong effects of flight distance (negative 
relationship between distance and susceptibility), breeding 
sociality (solitary species having greater susceptibility 
than colonial species), food (species eating live prey 
being less susceptible than herbivores), prevalence of 
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Tab. 2 Best-fit models of the relationship between susceptibility to mortality due to traffic and (models A-B) habitat, log- 
transformed nestling period and log-transformed body mass, (models C-D) square-root arcsine-transformed 
prevalence of Leucocytozoon and Plasmodium and log-transformed body mass, (models E-F) log-transformed mass 
of bursa of Fabricius and log-transformed body mass 





Variable Sum of squares df F P Slope (SE) Effect size 
Model A (Species): 
Habitat 535 1 9.14 0.0041 0.12 (0.04) 0.41 
Nestling period 6.03 1 10.30 0.0025 —1.20 (0.37) 0.43 
Body mass 8.00 1 13.67 0.0006 0.41 (0.11) 0.48 
Error 26.33 45 
Model B (Contrasts): 
Habitat 0.76 1 2.01 0.16 0.10 (0.07) 0.20 
Nestling period 1.84 1 4.89 0.032 —1.40 (0.63) 0.31 
Body mass 0.31 1 0.82 0.37 —0.22 (0.24) 0.13 
Error 17.38 46 
Model C (Species): 
Leucocyto-zoon 11.00 1 16.92 0.0002 0.72 (0.18) 0.51 
Plasmodium 6.66 1 10.25 0.0025 0.84 (0.26) 0.42 
Error 30.55 47 
Model D (Contrasts): 
Leucocyto-zoon 3.22 1 10.27 0.0027 0.84 (0.26) 0.46 
Plasmodium 3.85 1 12.25 0.0012 1.27 (0.36) 0.49 
Body mass 2.88 1 9.16 0.0044 —0.76 (0.25) 0.44 
Error 12.24 39 
Model E (Species): 
Bursa of Fabricius 13.93 dl 18.61 < 0.0001 0.41 (0.10) 0.55 
Error 32.18 43 
Model F (Contrasts): 
Bursa of Fabricius 1.88 1 4.46 0.041 0.57 (0.27) 0.31 
Body mass 2.70 1 6.42 0.015 —0.92 (0.36) 0.36 
Error 17.66 42 


Models A, C and E were based on species-specific data, while models B, D, and E were based on statistically independent contrasts. The overall models had the 
following statistics: Model A: F = 12.03, df= 3,45, r= 0.45, P < 0.0001, Model B: F = 2.82, df= 3,46, r= 0.15, P = 0.049, Model C: F = 13.13, df= 2,47, r 
= 0.36, P < 0.0001, Model D: F = 6.13, df= 3,39, 7? = 0.32, P = 0.0016, Model E: F = 18.62, df= 1,43, 7? = 0.30, P < 0.0001, and Model F: F = 3.28, df= 2,42, 
17? =0.14, P = 0.047. Effect size was estimated as the Pearson product moment correlation coefficient. 
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Fig. 3 Susceptibility to traffic death in different species of birds in relation to (A) breeding sociality, (B) relative prevalence 
of Plasmodium (after adjusting for prevalence of Leucocytozoon and body mass), (C) mass of the Bursa of Fabricius (g), 
and (D) flight distance (after adjusting for breeding sociality and mass of bursa of Fabricius) 


The lines are the linear regression lines weighted by sample size. 
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Tab. 3 Best-fit models of the relationship between susceptibility to mortality due to traffic and habitat, log-transformed 
nestling period, square-root arcsine-transformed prevalence of Leucocytozoon and Plasmodium, log-transformed 
mass of bursa of Fabricius and log-transformed body mass 








Variable Sum of squares df F P: Slope (SE) Effect size 
Model A (Species): 
Breeding sociality 9.83 1 19.77 < 0.0001 —0.38 (0.08) 0.59 
Flight distance 8.28 1 16.66 0.0002 —0.76 (0.19) 0.56 
Bursa of Fabricius 13.90 1 27.96 < 0.0001 0.49 (0.09) 0.66 
Error 18.40 37 
Model B 
(Contrasts) : 
Flight distance 2.50 10.43 0.0029 —0.93 (0.29) 0.50 
Breeding sociality 3.12 13.00 0.0010 —0.52 (0.14) 0.54 
Food 1.29 5.36 0.027 —0.33 (0.14) 0.38 
Plasmodium 2.94 1 12.24 0.0014 1.21 (0.35) 0.53 
Bursa of Fabricius 3.42 14.25 0.0007 0.89 (0.24) 0.55 
Body mass 2.67 11.12 0.0022 —1.15 (0.34) 0.51 
Error 7.68 32 


Model A was based on species-specific data, while model B was based on statistically independent contrasts. The overall models had the following statistics: 
Model A: F = 18.08, df= 3,37, r= 0.59, P < 0.0001, and Model B: F = 7.57, df=7,32, 1? =0.19, P < 0.0001. Effect size was estimated as the Pearson product 


moment correlation coefficient. 


Plasmodium (species with higher prevalence being more 
susceptible), bursa of Fabricius (species with larger bursa 
being more susceptible), and body mass (larger species 
being less susceptible) (Tab. 3 Model B). 


3 Discussion 


We provide evidence of avian mortality due to 
traffic directly reflecting the abundance of different 
species, but also phenotypic characteristics relating to 
flight behavior, visual acuity and health status. The null 
hypothesis investigated was that the abundance of traffic 
casualties simply reflected the abundance of different 
species, but also their behavior, with species with many 
individuals sitting on the road or flying close to the 
ground experiencing a high risk of mortality due to 
traffic. Abundance and to a small extent the proportion 
of individuals that were sitting on the road explained 
variation in traffic casualties. In addition, risk-taking 
behavior when approached by a human being estimated 
in terms of flight distance, breeding sociality (solitary 
species run higher risks), prevalence of the blood 
parasites of the genus Plasmodium (species with high 
prevalence run higher risks), and the size of the bursa of 
Fabricius (species with large immune defense organs run 
higher risks) all predicted susceptibility to traffic. We 
will briefly discuss these findings. 


Mortality among birds due to traffic may constitute 
a non-negligible fraction of overall mortality, and 
Svensson (1998) suggested that this factor is seriously 
under-estimated. Hansen (1982) estimated that 3.27 
million birds died on roads in Denmark in 1980-1981. 
Given that the total population of birds, including 
fledglings, consists of 16.29 million pairs (based on the 
mean estimates reported by Grell (1998)), and assuming 
that mean clutch size is 4.5 eggs (the mean for 512 
species breeding in the Western Palearctic (Cramp & 
Perrins, 1977—1994), the mean number of clutches per 
year being 1.6 (the mean for 512 species breeding in the 
Western Palearctic (Cramp & Perrins, 1977—1994)), and 
mean nest predation rate being 24% for 72 species 
1977-1994)), then the total 
population can be estimated as (16.29 million pairs X 2 


(Cramp & Perrins, 


individuals per pair) + (16.29 million pairs 4.5 eggs 
per clutch X 1.6 clutches per year X 0.76) = 121.72 
million individuals. Therefore, mortality due to traffic 
can be estimated as (3.27 million / 121.72 million), or 
2.7%. Given that mean adult mortality rate for 172 bird 
species from the Western Palearctic (Cramp & Perrins, 
1977-1994) was 38%, mortality due to traffic accounted 
for 2.7% / 38% = 7.1% of overall mortality. This 
proportion is bound to be an over-estimate because 
juvenile mortality is likely to be greater than 38%. 
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However, the conclusion remains that avian mortality 
due to traffic constitutes a non-negligible fraction of 
overall mortality in the order of 5%—-10% of overall 
mortality. Whether this has any effect on population size 
will depend on whether mortality caused by traffic acts 
in a density-dependent manner. To the best of our 
knowledge this provides the first estimate of the relative 
importance of traffic casualties for overall mortality in 
any group of animals. 

The null hypothesis for traffic casualties is that risk 
of mortality should reflect the abundance of different 
species and their flight behavior. Abundance explained 
42% of the variance in traffic mortality among species. 
Surprisingly, we found little evidence of flight behavior 
affecting the abundance of traffic casualties. There was a 
weak positive relationship between the proportion of 
individuals of different species being recorded on the 
road and road-kills, while the proportion of individuals 
flying at a height of less than 2 m did not predict traffic 
casualties. This leaves 60% of the variance to be 
explained by other factors such as sampling effort, 
ecological factors and other sources of variation. 

Blumstein (2006) suggested in his hypothesis about 
the evolutionary ecology of fear that flight distance 
should provide a reliable metric reflecting the risk that 
individuals are willing to take when trading foraging 
against flight. Studies of flight distance in birds have 
indicated that urban birds take greater risks than rural 
conspecifics, as reflected by shorter flight distances 
(Moller, 2008a). A number of factors hypothesized to 
provide increased visual acuity such as sociality and 
reliance on live prey may also allow individuals to 
escape when approached by a car. We found clear 
evidence suggesting that flight distance and breeding 
sociality both predicted susceptibility to traffic, and that 
was also the case in analyses of contrasts. 

Risk of mortality due to traffic was elevated in bird 
species with high prevalence of blood parasites of the 
genus Plasmodium. While we have no information 
whether individuals killed by cars were infected by blood 
parasites, we can conclude that species that had high 
prevalence of Plasmodium suffered disproportionately 
from traffic mortality. We have no a priori explanation 
why there were no effects for the blood parasite genera 
Haemoproteus and Trypanosoma, although previous 
studies of blood parasites and predation also showed 
similar variation among taxa (Moller & Nielsen, 2007). 
We hypothesize that this effect was due to the impact of 
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malaria as a cause of dramatic increases in the abundance 
of damaged erythrocytes, thereby severely reducing the 
efficiency of metabolism, but also restricting blood flow 
to the brain and the spleen (Atkinson & Van Riper, 1991; 
Chen et al, 2001). In addition blood parasites cause 
immune responses, fever and induction of heat-shock 
proteins (Wakelin, 1996; Merino et al, 1998; Chen et al, 
2001), and these effects may affect the ability of an 
individual to escape when approached by a car. Such 
effects will be particularly strong for intra-erythrocytic 
blood parasites like Plasmodium that is the cause of 
severe anemia. Bird species with blood parasite 
infections have shorter flight distances than uninfected 
species (Meller, 2008b), suggesting that species with 
high level infections may take greater risks when 
encountering potential danger. The increase in mortality 
risk by traffic linked to blood parasite infections selects 
against blood parasite infections. The findings reported 
here also have implications for parasitological studies of 
blood parasites. Given that traffic-induced mortality 
disproportionately biases estimates of prevalence of 
blood parasites downwards, any epidemiological 
inference based on prevalence in populations of hosts is 
likely to be biased likewise. 

Given the relationship between blood parasite 
infection and susceptibility to traffic, we should also 
expect immunity to explain interspecific variation in risk 
of being killed by a car. Indeed, we found a positive 
relationship between susceptibility to car accidents and 
size of the bursa of Fabricius. The bursa is an important 
immune defense organ involved in production of 
antibodies and responsible for differentiation of the 
entire repertoire of B-cells in young birds (Glick, 1983, 
1994; Toivanen & Toivanen, 1987). This effect was 
independent of body mass. 

Susceptibility traffic 
increased with increasing body mass, when analyzing 


to mortality caused by 


species-specific data and contrasts. Large bird species 
need longer distances for take-off, and climb rates 
decrease with body mass (Hedenström & Alerstam, 
1992), and species with large body mass have long flight 
distances (Meller, 2008a,b). That may partly explain 
why species with large body mass run a higher risk of 
mortality due to traffic than small species, although other 
factors relating to body size may explain this variation. 
The best-fit model showed a negative relationship with 
body mass still being present after other predictor 
variables had been included. 
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We also investigated susceptibility to traffic in 
relation to habitat, flight behavior, duration of development 
and several other factors, but found little evidence of 
these factors explaining interspecific variation once other 
factors were included in the phylogenetic analyses. 

In conclusion, we have found that risk of mortality 
in birds due to traffic increased with abundance. Once 
abundance had been taken into account, solitary species 
with short flight distances having a high prevalence of 
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Appendix 1 Frequency of road-kills of 50 species of birds, abundance of these species during standard point count censuses, proportion of 
individuals recorded on the road, proportion of individuals flying below 2 m, duration of nestling period (days), habitat (0 — grassland, 1 — 
bushes, 2—trees), breeding sociality (0 — solitary, 1 — colonial), mean flight distance (m), body mass (g), prevalence of Leucocytozoon and 


Plasmodium, and mass of bursa of Fabricius (g) 
































Species ofroads AP" "individuals, individuals, period F Breeding gisgance BOY Se renaney-PPEvsenee Of para of 
kills dance onroad flying below2m (d) bitat sociality (m) ass (g) tozoon Plasmodium Fabricius (g) 

Acrocephalus scirpaceus 1 0 0.00 1.00 0.9 0 0 6.73 11.80 0.01 0.02 0.016 
Actitis hypoleucos 1 0 1.00 1.00 27.0 0 0 17.00 47.75 = = = 

Alauda arvensis 0 25 0.90 1.00 90 0 0 13.92 36.40 0.00 0.00 0.050 
Ardea cinerea 0 1 0.00 0.00 50.0 0 1 47.36 1433.00 = = 1.663 
Buteo buteo 0 7 0.29 0.40 52.5 1 0 51.07 806.50 0.05 0.00 0.663 
Carduelis cannabina 6 5 0.50 0.00 3.5 1 0 10.87 18.95 0.00 0.01 0.045 
C. chloris 4 15 0.70 0.00 44 2 0 6.82 27.65 0.00 0.07 0.024 
Columba livia 3 0 1.00 0.00 36.0 0 1 3.50 261.00 0.00 0.00 0.197 
C. palumbus 7 0.38 0.60 33.5 2 0 12.14 494.50 0.13 0.14 0.049 
Corvus corone corone 1 12 0.84 0.14 32.2 2 0 74.00 544.50 0.33 0.01 0.656 
C. c. cornix 2 0 0.84 0.14 322 2 0 27.19 544.50 0.33 0.01 0.656 
C. frugilegus 2 3 0.48 0.05 330 2 1 34.17 453.50 0.00 0.00 0.764 
Cuculus canorus 0 1 0.00 0.67 9.0 1 0 25.57 120.50 0.00 0.00 0.068 
Delichon urbica 1 5 0.00 0.60 270 0 1 6.34 9.55 0.00 0.00 = 

Dendrocopos major 0 1 0.00 0.00 22.0 2 0 11.64 89.65 0.00 0.00 0.072 
Emberiza citrinella 8 12 0.83 0.33 2.0 1 0 12.51 26.75 0.05 0.03 0.053 
Erithacus rubecula 4 2 1.00 1.00 34 2 0 4.99 6.35 0.04 0.01 0.027 
Fringilla coelebs 14 23 0.84 0.88 3.9 2 0 8.30 24.20 0.14 0.00 0.036 
Hippolais icterina 2 6 0.00 0.00 3:5: 2 0 8.05 3.30 0.00 0.07 0.018 
Hirundo rustica 8 39 0.03 0.50 95 0 1 10.19 9.10 0.01 0.00 0.032 
Lanius collurio 1 0 0.00 0.00 4.5 1 0 6.10 30.70 0.02 0.03 0.045 
Larus argentatus 0 1 0.00 1.00 37:5 0 1 = 895.00 = # 0.862 
Larus canus 1 7 0.00 1.00 35.0 0 1 F 386.50 a z = 

L. fuscus 0 1 0.00 1.00 350 0 1 25.00 817.50 = = = 

L. ridibundus 0 4 0.00 1.00 350 0 1 R 280.50 = = 0.514 
Luscinia luscinia 1 0 0.00 1.00 9.6 1 0 15.89 25.00 = a = 

Motacilla alba 2 4 0.76 0.74 37 0 0 11.65 20.75 0.39 0.00 0.054 
Parus caeruleus 4 0 0.00 0.80 9.0 2 0 5.59 11.75 0.16 0.01 0.027 
P. major 4 22 0.07 0.57 89 2 0 5.32 18.50 0.03 0.04 0.030 
P. palustris 2 4 0.00 0.70 85 2 0 6.88 11.90 0.00 0.00 0.027 
Passer domesticus 199 261 0.83 0.85 4.1 0 1 3.39 30.35 0.00 0.07 0.043 
P. montanus 49 65 0.80 0.90 7.5 1 1 4.91 21.70 0.00 0.13 0.028 
Phasianus colchicus 3 5 0.67 .00 75.0 1 0 62.20 1389.10 0.15 0.50 0.157 
Phoenicurus ochruros 0 0.00 .00 5.5 0 0 6.78 16.00 0.00 0.00 0.006 
Phylloscopus collybita 11 0.00 0.80 50 2 0 6.71 7.70 0.04 0.00 0.011 
P. trochilus 0 2 0.00 0.20 323 2 0 8.06 9.35 0.04 0.01 0.009 
Pica pica 0 0.67 .00 27.0 2 0 14.16 228.00 0.45 0.00 0.470 
Prunella modularis 3 0 0.00 0.80 1:5 -2 0 5.09 18.95 0.06 0.00 0.045 
Pyrrhula pyrrhula 0 1 0.09 0.00 50 2 0 5.90 31.05 0.05 0.00 0.019 
Streptopelia decaocto 1 0.83 .00 7.0 2 0 5.36 201.50 0.00 0.00 0.045 
Strix aluco 0 0.00 0.00 34.5 2 0 = 491.00 1.00 0.00 0.442 
Sturnus vulgaris 1 24 0.12 0.10 210 2 1 10.77 80.50 0.0 0.00 0.118 
Sylvia atricapilla 3 0 0.00 0.50 10 2 0 5.73 18.85 0.0 0.00 0.044 
S. borin 2 1 0.00 0.50 1:2. -2 0 6.18 19.05 0.04 0.01 0.023 
Sylvia communis 5 11 0.00 0.80 1.0 1 0 7.95 14.50 0.0 0.00 0.016 
S. curruca 7 0.20 0.80 1.0 1 0 4.70 12.40 0.00 0.03 0.029 
Troglodytes troglodytes 4 0.10 1.00 73 2 0 5.51 8.90 0.0 0.00 0.048 
Turdus merula 95 35 0.36 0.83 3.6 2 0 743 95.85 0.1 0.13 0.237 
T. philomelos 2 3 0.40 0.80 32 2 0 8.96 70.50 0.30 0.02 0.153 
T. torquatus 1 0 0.20 0.80 50 2 0 S 117.00 0.44 0.00 0.132 





